Improving jet distributions with effective field theory 



\0. 

o: 

>: 
o. 

:z;; 

oo ■ 

(N ■ 

>: 

in . 

o. 

Oh: 

I 

Oh: 



X 



Christian W. Bauecl and Matthew D. Schwart4!l 
Ernest Orlando Lawrence Berkeley National Laboratory and University of California, Berkeley, CA 94720 

(Dated: February 2, 2008) 

We obtain perturbative expressions for jet distributions using soft-coUinear effective theory 
(SCET). By matching SCET onto QCD at high energy, tree level matrix elements and higher 
order virtual corrections can be reproduced in SCET. The resulting operators are then evolved to 
lower scales, with additional operators being populated by required threshold matchings in the effec- 
tive theory. We show that the renormalization group evolution and threshold matchings reproduce 
the Sudakov factors and splitting functions of QCD, and that the effective theory naturally com- 
bines QCD matrix elements and parton showers. The effective theory calculation is systematically 
improvable and any higher order perturbative effects can be included by a well-defined procedure. 



Calculating the distribution of jets in collider experi- 
ments is essential for understanding the standard model 
and in looking for signals of new physics. Unfortunately, 
improving the accuracy of a calculation is not simply 
a matter of computing higher loop QCD diagrams by 
brute force. These calculations produce large logarithms, 
which must be resummed. Moreover, fixed order calcu- 
lations are only feasible for a modest number of partons. 
For comparison with experimental data, however, it is 
desirable to obtain theoretical predictions for processes 
with a much larger number of final-state particles. The 
traditional method is to produce the required distribu- 
tions using splitting functions, which are derived in the 
limit of small transverse momentum. In this limit, dis- 
tributions factorize into lower order processes, multiplied 
by the classical probabilities for particles to branch. The 
probability for no branching is given by a Sudakov factor, 
which is simply an integral over the splitting functions, 
and sums the leading logarithmic (LL) terms. Several 
programs use this method of parton showers to generate 
distributions with an arbitrary number of final states and 
are valid at LL [l|, \^ . Since the underlying assumption 
in parton showers is small transverse momentum, addi- 
tional information from QCD matrix elements is required 
to properly describe jets which are widely separated. 

The issue of consistently matching fixed order QCD 
calculations with parton shower evolution has been the 
source of active research over the past several years [3j 
y, |5|. One must avoid double counting between emis- 
sions described by QCD matrix elements and those con- 
tained in the splitting functions, and cancel the infrared 
divergences between real and virtual QCD diagrams. 
Most programs include some hard matrix element correc- 
tions [6| , however there is no systematic way to improve 
their accuracy. 

Fundamentally, the problem of jet distributions is one 
of scale separation. These scales include the hard scale 
Q of the underlying interaction, and the scales p^ for 
the relative transverse momenta of additional partons. 
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Thus, it seems natural to reformulate the problem in the 
language of effective field theory, where scales are natu- 
rally separated. The appropriate effective theory for this 
problem is the soft-collinear effective theory (SCET) [3]. 
As we will show, the information of fixed order QCD is 
contained in SCET by a matching calculation at the hard 
scale Q, while the summation of the large logarithms is 
achieved by renormalization group (RG) evolution within 
SCET. Furthermore, in contrast to other approaches, the 
effective theory allows us to systematically incorporate 
higher order effects by straightforward and well-defined 
calculations. 

The matching at the hard scale involves determining 
the Wilson coefficients of SCET operators such that ma- 
trix elements in the effective theory reproduce the full 
theory up to a fixed order in as ■ Because all the infrared 
physics of the full theory is reproduced in SCET, the 
matching coefficients will be finite at any order in ag . To 
reproduce QCD with TV partons requires operators with 
up to N fields in the effective theory. To calculate this 
matching, of course, requires knowledge of QCD matrix 
elements with up to N partons in the final state. 

After the matching, full QCD is no longer needed. The 
next step is to evolve the operators in SCET to lower 
scales using the RG. This sums the large logarithms and, 
as we will show, it reproduces the Sudakov factors at LL. 
Lowering the renormalization scale in SCET corresponds 
to lowering the value of pt accessible to fields in the 
effective theory. As the scale gets lowered below the pt 
of one of the partons in the final state, the emission of 
this parton is no longer described by the interactions of 
SCET. To ensure that the theory still describes the same 
physics, a threshold matching is required which gives rise 
to new operators with additional coUinear fields. This 
matching involves SCET matrix elements for emission, 
which give splitting functions in the coUinear limit. 

For more details, we need to understand the basic con- 
struction of SCET. We want to describe the long distance 
behavior of QCD matrix elements with several energetic 
partons in well-separated directions. To achieve this, 
SCET requires coUinear fields for each of these directions, 
as well as soft fields they can interact with. CoUinear 
fields £,n are labeled by a light-like direction n^ and have 
energy much larger than the momentum transverse to n'^. 



E ^ p-^. CoUinear fields in SCET appear in operators 
as jets Xn which are fields wrapped in Wilson lines 



Xn = Wn£,n ■ 



(1) 



The Wilson lines ensure gauge invariance of the theory. 
For example, a two jet operator is 



O2 = XnFx,-, = LW^rWnXn 



(2) 



where F denotes some Dirac structure. To correctly de- 
scribe full QCD with up to N jets in different directions 
requires operators in SCET with up to N collinear fields, 
each with their own label n^. 

Since the sum of a number of collinear momenta is still 
collinear, there are interactions in SCET among parti- 
cles collinear to each direction nf . Particles collinear to 
different directions do not interact directly, but can ex- 
change soft gluons. The RG evolution in the theory, from 
higher to lower p^ , resolves the transverse momentum of 
jets. It physically disassembles jets into their partonic 
constituents. 

In this paper, we will study as an example the canon- 
ical process e+e~ —^ hadrons (for more details, see [iQi]). 
with center-of-mass energy denoted by Q. As mentioned 
above, the full matrix elements are matched onto SCET 
at the hard scale Q, while the RG evolution in the effec- 
tive theory sums the logarithms. To begin, we only in- 
clude the QCD matrix element with 2 partons in the final 
state, which allows us to focus on running. In this case, 
additional partons arise from SCET emissions through 
the threshold matching. Thus, the results will will only 
be valid in the limit Q ^ pj, ^ prp 3> ... , where pj, 
denote the transverse momenta between particles in the 
final state. Note that this is the same limit in which 
parton showers are known to hold. Because the effective 
theory includes the RG evolution, it will correctly sum 
all the logarithms of ratios of these scales. 

The matching condition is 



(QCD|99-)q = [C2(Ai)(02|gg) 



t^=Q- 



where {Q(JD\qq) denotes the fixed order matrix element 
of the current qTq. For particular final state momenta 
g^, we define ^n\q) = if n'^ is not aligned with q^^ . That 
is, the matrix element is non-zero only if n^ = q^ / Eq. 
Then the matching is trivial and 



C2(Q) = 1 



(4) 



Having performed this matching, we can focus on the 
resummation of the large logarithms. This is achieved by 
running the operator O2 below the scale Q. The Wilson 
coefficients satisfy the RG equation (ilogC„(//)/(ilog/i = 
7„(/i), which has the solution 



C„(/ii) _^ , . 

-r-f — 7 = n„(/i2,Mi) = cxp 



Ml /^ 



(5) 



n„ is an evolution kernel which is similar to the Sudakov 
factor in traditional parton showers. The anomalous di- 
mension 72 can be calculated by simple one-loop dia- 
grams in SCET [sj, giving 



I2W = Cf 



log- 



-M 



(6) 



where Cf = 4/3. Note that it is only the logarithmic 
term in the anomalous dimension which gives rise to the 
LL resummation. The constant term 3/2 gives rise to 
a subset of the subleading logarithms and can thus be 
omitted if only LL accuracy is desired. To obtain the 
complete NLL result, the coefficient of the log n/Q term 
is required at two loops. Higher orders in the resumma- 
tion can straightforwardly be included by calculating the 
anomalous dimensions at even higher loop order. 

At a scale px a third parton can be resolved in one 
of the jets. Below that scale the matrix element of O2 
with the three partons vanishes and a threshold match- 

f21 

ing onto a new operator Og has to be performed (the 
superscript indicates that this operator arises through 
threshold matching). The matching condition is 



M=PT+e 



cf)(A.)(0f|gg-5) 



[C2{fi){02\qqg)] 

(7) 
Emissions in SCET can be collinear or soft, and the 
collinear gluon can either be radiated from the collinear 
quark (anti-quark), £,ni£,n) , or it can be emitted by one of 
the Wilson lines, W^(T^fi)- There are some ambiguities 
in calculating the matrix elements in SCET, which will 
only affect higher order contributions if incorporated con- 
sistently. Here, we will take the additional gluon in the 
final state to be collinear, since soft gluons have much 
smaller p-^ than collinear gluons. We also use the fact 
that the emission of a collinear gluon from the collinear 
Wilson line will give rise to longitudinally polarized glu- 
ons (up to higher orders in pr/Q)- Thus, to calculate the 
matrix element of O2 in Eq. ([7]) we include only the emis- 
sion from the collinear fermions and keep only transverse 
(3) gluon polarizations e^. With these conventions Eq. ^ 



is satisfied by C3 '{pt) — C2{pt) and 



{o'i^\qqg) = gsXn, 



^ 



±^ql 



F%/J 



T.q-{Pq+Pg) nq-{pq+pg) 



Xn, .(8) 



Note that we choose to integrate out emissions from 
both the quark and the anti-quark at the same scale pr- 
In general, the transverse momentum of the gluon with 
respect to the quark can differ from the one with respect 
to the anti-quark, and one could choose to integrate out 
each emission at these different scales. However, the dif- 
ference between these scales is subleading in pr/Q, and 
thus beyond the order we are working. 

This procedure of matching and running is continued 
for each value of p^ . A closed form expression can be 
obtained for the anomalous dimension of a general oper- 
ator On which contains Uq quark and Hg gluon fields. At 



ing order 






[2 2 J 



(9) 



and the solution to the RG equation is given by Eq. ^. 
At leading log accuracy the evolution kernels for the Wil- 
son coefScients are directly related to the Sudakov factors 



for quarks and 
parton showers 



luons, Ag and Ag, used in traditional 



n„(Ai2,Ml) = A^'/2(^2,Mi)A^''/'(M2,m). (10) 

The threshold matching can be continued for addi- 
tional emissions. In general, the rule for quark splitting 
in SCET can be written as 



Xr. 



i H' 



"2n-{pq+pg) 



Xr, 



(11) 



where the +(— ) is for fermions (antiferniions). A sim- 
ilar equation can be obtained for giuons splitting into 

(2) 
two giuons or two fermions. New operators Om are thus 

(2) 
obtained by replacing each coUinear field in 0,„_i with 

Eq. pT|) and summing over all contributions. Note that 

(2) 
the operators Om are sums over many different terms, 

each giving different splitting histories, just as the oper- 

ator O3 in Eq. ([8]) is the sum of splittings of the quark 

and the anti-quark. 

This sequence of threshold matching and running is 

repeated until a low scale /i = /io is reached, where non- 

perturbative physics becomes important. At that scale 

the amplitude of e^e~ — > hadrons is 



(SCET), 



J2C^^Hf^o){0^^^),o, 



(12) 



where the evolution of the operators between the various 
threshold matching scales is captured in Cm (/-to) 

c(^)(/io) = C2(Q)n2(Q,p«) • • • n„_i(p(r-^,Vo) . (13) 

To compute differential cross sections, we need to 
square this amplitude and sum over final spins of the 

fermions and transverse polarizations of the giuons. Be- 

(2) 
cause the operator Om is obtained from the operator O2 

by iterating Eq. (jlip. we can simplify the form of \Om p. 
To leading order in pt/Q, the result is 



E 

spins, pols 



0(2) 



spins, pols 



(2) 



Pt 



where P{z) is a splitting function, and the variable z de- 
notes the ratio of the energy of the emitted particle to the 
energy of the particle it was emitted from. Each emission 
of an additional particle gives rise to a splitting function 
after the matrix element is squared. Using Eq. P^ to- 
gether with Eq. (fTU]) , the cross section obtained from the 



amplitude (SCET) given in Eq. p2| reduces to the cross 
section of e+e^ -^ qq scattering, multiplied by products 
of splitting functions and Sudakov factors. Thus, SCET 
is equivalent to the traditional parton shower at leading 
order. 

Now that we have seen how SCET sums the large loga- 
rithms and showers additional coUinear particles, we can 
go back and incorporate the QCD matrix element with 
an additional gluon in the final state. Thus, we should 
include O2 and an operator O3 in the matching at /z = Q- 
This allows us to describe the differential rate with three 
well-separated jets, for any value ofprp, and will improve 
our previous results such that the resulting distributions 
are valid in the limit Q,prp ^ Pj, ^ pj, . . . . In order 
to obtain the correct normalization of the total cross sec- 
tion at 0{as), the one loop matching onto the operator 
O2 is required at the hard scale as well. 

To match onto O3 we calculate matrix elements in 
QCD and SCET with three partons in the final state 



{qCB\qqg) = C2{02\qqg) + CsiO^lqqg) . 



(15) 



The operator O2 and its Wilson coefficient C2 at order 
a^ have already been determined from the matching with 
two partons. With our convention for emissions 



{Oz\qqg) = gsXn 



/^^r 



pj^^j 



iq-iPq+Pg) iq-{Pq+Pg) 



Xn,(-16) 



and then Eq. p^ is satisfied for C3(Q) ~ 1. This matrix 
element vanishes in the limit that the gluon is soft or is 
coUinear to either the quark or the anti-quark. This is 
again due to the fact that SCET reproduces the IR of 
QCD. 

To obtain C2 to order Us we again use Eq. ([3]), but now 
calculate the matrix elements in QCD and SCET at one- 
loop. Since the IR divergences of QCD are reproduced 
in SCET, the Wilson coefficient is finite H 



C2(Q) = 1 



O-sCf 

Att 



6 



Stti 



(17) 



Even higher order corrections to these Wilson coefficients 
can be included by performing the matching at higher 
loop order. 

After the matching at fj, = Q, the operators are run 

down to a low scale. The running of O3 is identical to 

(2) 
the running of Og , since running in SCET depends only 

on the field content of an operator, not on its tensor 

structure. We run both O2 and O3 down to the scale 

(2) 
fi ^ Pt, at which point we match O2 onto Cg as before. 

At fi ^ Pt 



(SCET(3))p^ = C2 n2 (of) + C3 Ha {O: 



(18) 



After that, the sequence of running and matching con- 
tinues as before, but now with two series of operators: 

(2) 

Om , populated by branchings descended from O2, and 
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FIG. 1: Thrust distribution in 3-jet events, calculated as de- 
scribed in the text. Note how the SCET result reduces to 
QCD at small thrust (large pr) and to the parton shower at 
large thrust (small pr)- 



(3) 

Om , which are descended from O3. 
the scale no is 



The final result at 



(scet(3))^„ =n3(p^'\p^'')---n™-i(p^""'\Mo) x (i9) 
C2(Q)n2(Q,p«)(o(2)^ + C3(g)n3(Q,p^'^)(oi?) 

Note that this is very similar to what we had in Eqs. (fll 
and (fT3|) . But it is correct to 0{as) and for values of pt 
that satisfy Q , plj,' :$> plj.' :$> plj,' . . . . It reproduces the 
three jet differential cross section at NLO, but also sums 
all leading logarithms. 

The SCET result interpolates between fixed order 
QCD and a parton shower. To demonstrate this, it is 
enough construct observables based on 3-jet final states. 
Since we only need information about the first branch- 
ing we should evaluate the matrix elements at a scale 
of order pT, and the result is given by Eq. P^ . It 
is not hard to see that Eq. P^ has the right limits. 
For pt <S^ Q, SCET is an excellent approximation to 
QCD and the matrix element of O3, which represents 
the difference between QCD and SCET, vanishes. Also, 
1(03 \qqg)\'^ reduces to a splitting function in this limit, 



as indicated in Eq. ((14)) . and 112 reproduces the LL Su- 
dakov factor, as shown in Eq. pU)) . Thus, for pr ^ Q, 
SCET reduces to the product of a Sudakov factor and a 
splitting function, as in a parton shower. In contrast, for 
PT '- Q, SCET(3) is very close to QCD, since SCET has 
not evolved far from the scale where we matched it to 
QCD exactly. Given the results in these two limits, the 
effective theory will smoothly interpolate between QCD 
and the parton shower. To see this, in Fig. [T] we show the 
differential thrust distribution for Q = 91 GeV. For three 
partons, T = max{Eq,Eq,Eg}/2Q. The SCET result is 
given by squaring (jlSp and integrating over the remain- 
ing phase space, while for QCD we use the leading tree 
level differential cross section. For the parton shower, we 
use the splitting functions l/{pq -t-pg)^P(z), multiplied 
by the LL Sudakov factor. We have chosen the value of 
the scale ^ ^ px such that the maximum value of the 
Wilson coefficients and Sudakov factors is unity. 

We have shown that SCET reproduces QCD at NLO, 
including full parton showering. There are many ways 
the results presented here can be improved: 1) Include 
more operators at fj, = Q. 2) Include the matching at 
fi — Q to higher order in a^. 3) Include the running of 
the operators at subleading order. 4) Include the thresh- 
old matching in SCET at higher order. 5) Include power 
corrections from the SCET expansion. It is important 
to bear in mind that all these improvements are straight- 
forward calculations in the effective theory. For example, 
to include more operators at /i = Q simply requires ad- 
ditional tree level calculations in QCD and the effective 
theory, while to go beyond NLO only requires the calcu- 
lation of well-defined loop diagrams in QCD and SCET. 
No new formalism needs to be developed to achieve any of 
these results. Thus, SCET is a powerful tool for improv- 
ing the theoretical understanding of jet distributions. It 
provides a convenient systematically improvable frame- 
work for performing higher loop calculations, and should 
be straightforward to implement in an event generator, 
allowing direct comparison to data. 
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